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Abstract 

For  the  first  time,  the  dynamic  behavior  of  a  thermally  coupled  hydrogen  storage  and  fuel  cell  system  was  studied  using  experimentally 
validated  models  of  a  metal-hydride  hydrogen  storage  system  and  a  proton  exchange  membrane  (PEM)  fuel  cell  stack.  Two  specific  situations 
were  studied:  in  the  first  scenario,  the  fuel  cell  and  metal  hydride  bed  were  thermally  coupled;  in  the  second  scenario,  they  were  thermally 
uncoupled.  The  results  from  these  simulations  clearly  exposed  unique  and  subtle  behavior  associated  with  thermally  coupled  systems  that 
could  not  be  easily  gleaned  from  simulating  each  device  alone,  as  revealed  from  studying  the  uncoupled  case.  The  most  important  result 
exemplified  the  importance  of  providing  waste  (excess)  energy  from  the  PEM  fuel  cell  to  the  metal  hydride  bed  to  facilitate  the  removal  of 
hydrogen  during  a  discharge  process.  A  detailed  parametric  study  is  forthcoming. 

©  2004  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

The  use  of  hydrogen  as  an  alternative  to  carbonaceous  fu¬ 
els  is  of  increasing  interest.  Hydrogen  exhibits  very  attractive 
features  for  use  in  applications  such  as  hybrid  electric  vehi¬ 
cles  [1],  because  it  can  be  used  essentially  pollution-free  and 
produced  from  renewable  energy  resources,  thus  eliminat¬ 
ing  the  net  production  of  greenhouse  gases.  Recent  studies 
have  also  indicated  that  hydrogen  costs  are  reasonable,  which 
make  it  an  ideal  candidate  to  replace  fossil  fuels  as  an  energy 
source  [2]. 

Hydrogen  storage  has  been  the  subject  of  intensive  re¬ 
search  for  many  years.  There  are  currently  three  main  tech¬ 
nologies  for  storing  hydrogen  on  vehicles:  as  a  compressed 
gas,  as  a  liquefied  cryogenic  fluid,  and  as  a  solid  in  a  metal 
hydride  [3].  The  automotive  industry  has  given  considerable 
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attention  to  storing  hydrogen  in  metal  hydrides  because  it 
can  be  stored  reversibly  in  the  solid  state  at  relatively  low 
pressures  and  ambient  temperatures  [4-7].  The  storage  ves¬ 
sel  typically  contains  a  powdered  metal  (often  an  alloy)  that 
absorbs  hydrogen  exothermically  with  the  concomitant  re¬ 
lease  of  energy  as  the  tank  is  filled  with  hydrogen  under 
pressure.  The  metal  hydride  desorbs  the  hydrogen  endother¬ 
mically  by  reducing  the  pressure  and  supplying  some  energy. 
While  many  applications  take  advantage  of  these  properties, 
including  rechargeable  batteries  [8],  and  heating  or  cooling 
systems  [9],  this  paper  addresses  the  application  of  metal- 
hydride  hydrogen  storage  for  use  as  a  fuel  in  energy  conver¬ 
sion  systems  such  as  fuel  cells. 

Fig.  1  illustrates  a  typical  hydrogen-fueled  energy  system 
that  is  used  in  a  hybrid  electric  vehicle.  The  fuel  cell  stack 
together  with  the  boost  converter  comprises  the  power  gen¬ 
eration  unit.  The  battery  stores  energy  when  the  power  is 
sufficient  and  releases  energy  when  peak  power  is  needed. 
When  accelerating  the  vehicle,  both  the  fuel  cell  and  the 
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Nomenclature 

Acen  effective  area  of  each  cell  in  the  fuel  cell  stack 
(cm2) 

Aex  effective  area  of  the  fuel  cell  heat  exchanger 
(cm2) 

bj(j  =1,2)  Henry’s  law  constant  defined  in  Eq.  (4) 
(Pa- 1/2 ) 

bjo(j=  1,2)  temperature  independent  coefficient  of  the 
Henry’s  law  constant  (Pa- 1/2 ) 

Cpc  average  heat  capacity  of  the  fuel  cell  stack, 
(Jkg-1K-1) 

Cpn?  heat  capacity  of  H2  (J  mol  - 1  K- 1 ) 

Cps  heat  capacity  of  the  metal  hydride  ( J  kg- 1  K- 1 ) 

dj(j= 0, 1, 2)  coefficients  in  Eq.  (22) 

Do  coefficient  in  Eq.  (5) 

Dstack  voltage  of  the  fuel  cell  stack  (V) 

Eq  standard  potential  of  H2-O2  reaction  (V) 

/  defined  by  Eq.  (16) 

F  Faraday  constant 

h  simulation  step  time  (s) 

hex  overall  heat  transfer  coefficient  of  the  fuel  cell 
heat  exchanger  (J  m-2  K-1  s-1) 
hbed  overall  heat  transfer  coefficient  of  the  metal 
hydride  bed  heat  exchanger  (J  m-2  K-1  s-1) 
AH  isosteric  heat  of  adsorption  of  the  hydrogen  in 
the  metal  hydride  (J  mol-1) 

AH\  isosteric  heat  of  adsorption  below  the  switch 
pressure  (J  mol-1) 

AH2  isosteric  heat  of  adsorption  above  the  switch 
pressure  (J  mol-1) 

A TTS  coefficient  in  Eq.  (5)  (J  mol-1) 

A Hc  the  heat  of  combustion  of  hydrogen  (J  mol  - 1 ) 

I  the  current  from  the  fuel  cell  (A) 

kq  mass  transfer  coefficient  (s- 1 ) 

kp  defined  by  Eq.  (18) 

kcp  defined  by  Eq.  (19) 

L  length  of  the  hydride  bed  (m) 

rape  mass  of  the  fuel  cell  stack  (kg) 

Mmh  average  molecular  weight  of  metal  hydride 
(kg  mol-1) 

N  number  of  the  cells  in  the  fuel  cell  stack 
m  hydrogen  flow  rate  out  of  the  metal  hydride  bed 
(mol  s-1) 

nj(j= 0, 1, 2)  coefficients  in  Eq.  (22) 
h  hydrogen  flow  rate  (mol  s  - 1 ) 

P  pressure  inside  the  metal  hydride  bed  (Pa) 

P0  pressure  inside  the  metal  hydride  bed  at  the 
initial  state  (Pa) 

Pfc  pressure  of  the  hydrogen  in  the  fuel  cell  (Pa) 
Ps  equilibrium  pressure  of  the  hydrogen  in  the 
metal  hydride  bed  at  the  envelope  (Pa) 
q  hydrogen  concentration  in  the  metal  hydride 

(mol  mol-1) 


q  equilibrium  hydrogen  concentration  in  the 

metal  hydride  (mol  mol-1) 

gmax  maximum  equilibrium  hydrogen  concentration 

in  the  metal  hydride  (mol  mol-1) 
q$  equilibrium  hydrogen  concentration  in  the 
metal  hydride  at  the  initial  state  (mol  mol-1) 
qs  saturation  hydrogen  concentration  in  the 

metal  hydride  at  the  low  concentration  phase 
(mol  mol-1) 

Q  rate  of  energy  exchange  with  the  metal  hydride 
bed  (J  s-1  m-3) 

rj(j  =  0, 1, 2)  coefficients  in  Eq.  (22) 

R  ideal  gas  constant  (J  mol-1  K-1) 

Ri  radius  of  the  metal  hydride  bed  heat  exchanger 

tube  (m) 

R0  radius  of  the  outer  wall  of  the  metal  hydride 
bed  (m) 
t  time  (s) 

T  average  temperature  of  the  hydrogen  inside  the 

metal  hydride  bed  (K) 

Tamb  the  temperature  of  the  ambient  (K) 

To  average  temperature  of  the  hydrogen  inside  the 
metal  hydride  bed  at  the  initial  state  (K) 

Tm  temperature  of  the  heat  transfer  medium  for  the 

metal  hydride  bed  (K) 

Tpc  the  bulk  temperature  of  the  fuel  cell  stack  (K) 
Vb  volume  of  the  metal  hydride  bed  (m3) 

Wamb  thermal  power  transferred  from  the  fuel  cell  to 
the  ambient  environment  (W) 

We  electrical  power  released  from  the  fuel  cell  (W) 
Wth  thermal  power  generated  by  the  fuel  cell  (W) 

Greek  symbols 

s  metal  hydride  bed  porosity,  dimensionless 

ps  density  of  the  metal  hydride  (kg  m-3) 
pi  rate  of  pressure  change  inside  the  metal  hydride 

bed  (Pas-1) 

V2  rate  of  average  temperature  change  of  the  hy¬ 
drogen  in  the  metal  hydride  bed  (K  s-1) 

P3  rate  of  H2  concentration  change  inside  the 
metal  hydride  bed  (mol  mol-1  s-1) 


Fig.  1.  Block  diagram  of  a  typical  hydrogen-fueled  energy  system  that  is 
used  in  a  hybrid  electric  vehicle. 
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battery  provide  electrical  power  to  the  permanent  magnet 
(PM)  synchronous  motor  through  the  PWM  inverter.  When 
decelerating  the  vehicle,  the  motor  acts  as  a  generator  and  the 
generated  power  is  used  to  charge  the  battery.  The  hydrogen 
storage  system  provides  hydrogen  for  the  fuel  cell  stack  as 
necessary. 

When  hydrogen  is  released  from  the  metal  hydride  bed, 
the  endothermic  desorption  process  causes  its  temperature  to 
decrease  [10].  In  turn,  this  decrease  in  temperature  greatly  di¬ 
minishes  the  rate  at  which  hydrogen  can  be  released  from  the 
bed.  To  ameliorate  this  undesirable  effect,  some  of  the  excess 
heat  naturally  produced  from  the  fuel  cell  during  the  genera¬ 
tion  of  electricity  can  be  transferred  to  the  metal  hydride  bed 
to  increase  its  temperature.  In  other  words,  the  metal  hydride 
bed  can  be  thermally  coupled  to  the  fuel  cell  stack,  which  not 
only  facilitates  the  release  of  hydrogen  from  the  metal  hy¬ 
dride  bed,  but  it  also  helps  remove  excess  (waste)  heat  from 
the  fuel  cell  stack. 

The  objective  of  this  paper  is  to  study,  for  the  first  time, 
the  behavior  of  a  thermally  coupled  hydrogen  storage  and 
fuel  cell  system.  This  system  level  simulation  is  carried  out 
using  a  very  friendly  computational  environment  refereed 
to  as  the  virtual  test  bed  (VTB)  [11-14],  along  with  suit¬ 
able  models  that  describe  the  dynamic  performances  of  the 
metal  hydride  bed  and  fuel  cell  stack  with  enough  rigor  to 
capture  the  physics  but  not  overwhelm  the  system  simula¬ 
tor.  Although  many  detailed  models  of  a  metal  hydride  bed 
have  been  developed  to  study  the  dynamics  inside  the  bed 
[10,15-17],  only  that  developed  recently  by  Gadre  et  al.  [10] 
has  been  validated  against  extensive  experimental  data  and  is 
simple  enough  for  system  level  studies  [10].  For  similar  rea¬ 
sons,  a  relatively  simple  proton  exchange  membrane  (PEM) 
fuel  cell  model  is  chosen  here  for  these  system  studies  [18]. 

First,  the  inherently  nonlinear  model  equations  that  realis¬ 
tically  account  for  heat  and  mass  transfer  resistances  associ¬ 
ated  with  the  discharge  of  hydrogen  from  a  metal  hydride  bed 
[10]  are  reformulated  in  a  resistive  companion  (RC)  form  for 
efficient  solution  in  the  VTB  platform  [13].  Then  the  PEM 
fuel  cell  model  [18]  is  dynamically  coupled  to  the  RC  form 
of  the  metal  hydride  bed  model  in  the  VTB  environment. 
System-level  simulations  of  the  resulting  integrated  metal- 
hydride  hydrogen  storage  and  fuel  cell  system  are  then  car¬ 
ried  out  for  both  thermally  coupled  and  uncoupled  scenar¬ 
ios.  These  two  scenarios  expose  the  advantages  of  providing 
waste  or  excess  energy  from  the  PEM  fuel  cell  stack  to  the 
metal  hydride  bed  to  facilitate  the  removal  of  hydrogen. 


2.  Metal-hydride  hydrogen  system  model 


priate  assumptions  are  made  so  that  the  distributed  thermal 
and  mass  transfer  phenomena  that  occur  within  this  vessel 
are  described  by  lumped  approximations  that  still  allow  for  a 
realistic  but  less  complicated  description  of  this  complex  pro¬ 
cess.  First,  it  is  assumed  that  the  system  is  symmetric  around 
the  bed  axis  and  that  all  the  variables  vary  along  the  bed  ac¬ 
cording  to  plug  flow  conditions,  i.e.,  with  no  radial  gradients. 
Similarly,  axial  gradients  of  temperature  and  pressure  are  ne¬ 
glected  while  the  gas  molar  flux  in  the  bed  is  assumed  to  vary 
linearly  from  a  value  of  zero  at  the  closed  end  of  the  bed. 
Identical  temperatures  are  assumed  for  both  the  solid  and  gas 
phases.  The  external  surface  is  assumed  to  be  perfectly  insu¬ 
lated,  and  heat  transfer  between  the  bed  and  the  internal  heat 
exchanger  is  described  by  an  overall  heat  transfer  coefficient. 
The  temperature  inside  the  heat  exchanger  is  assumed  to  be 
constant  along  the  bed  and  throughout  the  hydrogen  desorp¬ 
tion  process.  The  mass  transfer  resistance  in  the  solid  phase 
is  described  by  a  linear  driving  force  (LDF)  approximation, 
using  a  mass  transfer  coefficient  that  is  assumed  to  be  inde¬ 
pendent  of  temperature.  The  desorption  process  starts  from  a 
fully  charged,  equilibrium  state  at  a  given  initial  temperature 
and  pressure.  Throughout  the  process,  it  is  assumed  that  the 
flow  rate  at  the  discharge  end  of  the  bed,  i.e.,  the  molar  flux,  is 
defined  either  as  a  constant  or  as  a  time  dependent  function. 

Based  on  these  assumptions  and  assuming  ideal  gas  be¬ 
havior,  the  one-dimension  mass  balance  is  governed  by 

1  dP  P  dT  h  l  —  £  Ps  dq 

- { - { - — - -  =  0 

RT  dt  RT 2  dt  n(Ro  —  R\)L  s  2Mmk  dt 

(1) 

where  P  is  the  pressure  inside  the  bed,  T  the  average  tempera¬ 
ture  of  the  hydrogen  in  the  metal  hydride  bed,  R  the  universal 
gas  constant,  h  the  hydrogen  flow  rate  demand,  Rx  the  radius 
of  the  heat  exchanger  tube,  R0  the  radius  of  the  bed,  L  the 
length  of  the  bed,  s  the  porosity  of  the  metal  hydride  bed,  ps 
the  density  of  the  metal  hydride,  Mmh  the  average  molecular 
weight  of  the  metal  hydride,  and  q  is  the  hydrogen  concentra¬ 
tion  in  the  metal  hydride.  The  corresponding  energy  balance, 
including  expansion/compression  effects,  but  ignoring  con¬ 
vection  effects  (as  it  contains  no  spatial  gradients),  is  written 
as 


Cpu2  + 


3  T  3  P  1  -  s  ps  A  H  3  q 
dt  dt  £  2Mmh  dt 


+ 


2/Zbed  Ri 

s{Rl  -  R() 


(T  -  Tm)  =  0 


The  metal  hydride  bed  under  study  consists  of  a  cylin¬ 
drical  container  of  radius  R0  filled  with  a  metal  powder  that 
chemically  absorbs  the  hydrogen  under  pressure.  This  con¬ 
tainer  also  has  an  internal  co-axial  heat  exchanger  of  radius  R[ 
that  supplies  the  energy  needed  to  release  the  hydrogen  during 
discharge  and  carries  away  excess  heat  during  charge.  Appro- 


where  Cp h7  is  the  heat  capacity  of  the  hydrogen,  Cps  the  heat 
capacity  of  the  metal  hydride,  AH  the  isosteric  heat  of  adsorp¬ 
tion  of  the  hydrogen  in  the  metal  hydride,  /zbed  the  overall  heat 
transfer  coefficient  of  the  metal  hydride  bed  heat  exchange, 
and  Tm  is  the  temperature  of  the  metal  hydride  bed  heat  ex¬ 
changer  exchange  medium.  The  intraparticle  mass  transfer 
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mechanism  is  based  on  the  following  LDF  expression: 


dq 

f  =  kq(q  -  q)  (3) 

ot 

where  q  is  the  concentration  of  hydrogen  in  the  metal  hydride 
in  equilibrium  with  the  gas  phase  at  P  and  T ,  and  kq  is  the 
mass  transfer  coefficient. 

The  thermodynamic  relationship  that  describes  the  equi¬ 
librium  loading  of  hydrogen  in  the  metal  hydride  as 
a  function  of  pressure  and  temperature,  the  so  called 
pressure-composition-temperature  (PCT)  relationship,  can 
be  adequately  represented  by  the  composite  Langmuir  model 
[10].  This  thermodynamic  model  consists  of  two  Langmuir- 
type  isotherms  connected  at  a  so-called  switch  point,  i.e., 


qs  (l  +  hVK) 

(l+biJPj  JPs’ 

(^max  q&)b2  (y/~P  a/Rs  ) 
qs  ^  7 - 7Z - \ — 

(l+b2(VP-  VK)J 


P  <  PS(T) 


p>  Ps(T) 


where  b\  =  Pi0exp  b2  =  /?20exp  A#i 

and  AH2  are  the  heats  of  adsorption  below  and  above  the 
switch  point,  and  qm ax,  b i0,  and  b2o  are  fitting  parameters.  Ps 
and  qs  are  the  switch  pressure  and  hydrogen  concentration  in 
the  metal  hydride  that  correspond  to  the  point  on  the  isotherm 
where  the  phase  transition  starts.  The  switch  variable  Ps  is 
related  to  T  according  to 


ln(Ps)  —  A)  + 


A//s 

RT 


where  D0  and  A Hs  are  additional  fitting  parameters.  Due 
to  the  composite  nature  of  the  Langmuir  isotherm  model, 
the  isosteric  heat  of  adsorption  of  the  hydrogen  in  the  metal 
hydride  is  evaluated  distinctively  depending  on  whether  the 
pressure  is  below  or  above  Ps ,  i.e., 


AH  = 


2AH{  =  A  Hs 


AH2  +  (AHx 


,  P  <  P*(T) 

,  P  >  PS(T )  (6) 


To  complete  the  description  of  the  model,  the  following 
initial  conditions  are  needed:  at  t-  0,  P  =  Po ,  T=Tq,  and 
q  =  qo  =  qo*- 

Eqs.  (1)— (6)  constitute  a  relatively  simple  mathematical 
model  that  realistically  accounts  for  heat  and  mass  transfer 
resistances  during  the  dynamic  discharge  of  hydrogen  from 
a  metal  hydride  bed.  It  should  be  noted  that  this  model  is 
a  simplified  version  of  a  complete  set  of  validated  equations 
that  describe  the  dynamics  of  the  hydrogen  discharge  process 
over  a  wide  range  of  hydrogen  flow  rates  [10].  The  main 
difference  is  that  axial  temperature  gradients  (i.e.,  convection 
effects)  are  ignored  in  this  simplified  version. 

This  resistive  companion  (RC)  method  was  applied  to 
the  system  of  equations  to  yield  a  form  of  the  model  that 
is  amenable  to  “natural  coupling,”  i.e.,  a  form  of  the  model 


that  can  be  inserted  into  a  network  without  regard  to  causal¬ 
ity  [13,14].  Natural  conservation  laws  are  automatically  en¬ 
forced  where  appropriate.  As  shown  above,  Eqs.  (1)— (3)  are 
nonlinear  differential  equations.  The  construction  of  the  RC 
model  for  the  metal  hydride  bed  begins  by  introducing  three 
intermediate  variables  to  replace  the  derivative  terms  in  the 
system  equations  by  defining 


Vl  = 

dt 

(7) 

dT 

v2  = 

dt 

(8) 

dq 

(9) 

V3  = 

dt 

where  v\,  v2,  and  v2  are  intermediate  variables.  Integrating 
Eq.  (7)  over  one  time  step,  i.e.,  from  time  instant  ( t  —  h)  to  t , 
allows  it  to  be  expressed  in  discrete  form  as, 

[1 nit)  + 1 nit  -h)]^  =  P(t)  -  Pit  -  h)  (10) 

where  t  is  the  current  simulation  time,  and  h  is  the  simulation 
time  step.  Rearranging  Eq.  (10)  yields 


0  =  P(t )  -  ^Vi it) 


Pit  -  h)  +  ^vi it  -  h ) 


(11) 


Similarly,  discretizing  Eqs.  (8)  and  (9)  and  rearranging 
gives  the  following  equations: 


0  =  7X0  -  \v2 (0  - 
h 

0  =  qit)  -  -v3(0  - 


h 

Tit  -  h)  +  -V2 it  -  h) 
h 

qit  -h)  +  -t>3 it  ~  h) 


(12) 


(13) 


Eqs.  (1 1)— (13)  are  in  the  standard  RC  form.  Based  on  the  def¬ 
initions  given  in  Eqs.  (7)-(9),  Eqs.  (l)-(3)  are  then  converted 
to  three  nonlinear  algebraic  equations  as  follows. 


kmv\ 

m  =  - 

RT 


km  P 

~rT1 


v2  +  kmkpV3 


(14) 


Q  —  kms 


C  ph2  +  kc 


v2  +  kms  v\  +  kms  kp  AH  v2 


(15) 


f  =  v 3  -  kq(q*  -q)  =  0 


(16) 


where  Q  is  the  rate  of  heat  exchange  between  the  metal  hy¬ 
dride  bed  and  its  heat  exchanger,  m  the  molar  flow  out  of  the 
hydride  bed,  V b  the  volume  of  the  metal  hydride  bed,  and  kp 
and  kcp  are  given  by  Eqs.  (17)— (21). 


Q  —  2/zbe(j7T  R[  L{Tm  T) 

k  _  1  ~ £  ^ 
p  s  2Mmh 

1  —  £ 

kcp  =  -  PsCAs 

c 


(17) 

(18) 

(19) 
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Fig.  2.  Model  predictions  of  the  experimental  pressure  histories  [10]  in  the 
metal  hydride  bed  for  hydrogen  flow  rate  demands  of  5  and  10  SLPM. 

km  =  jt(R20  -  R\)L  (20) 

m  =  —h  (21) 

The  RC  equations  are  written  in  matrix  form,  including  the 
conductance  matrix  and  history  vector.  In  this  RC  form,  the 
dynamic  model  of  the  metal  hydride  bed  is  easily  imple¬ 
mented  in  the  VTB,  with  seamless  integration  between  other 
components  and  devices  of  interest,  like  the  fuel  cell  stack. 

Figs.  2  and  3,  respectively,  compare  the  experimental  pres¬ 
sure  and  temperature  histories  [10]  with  those  obtained  from 
the  RC  form  of  the  metal  hydride  bed  model  for  hydrogen 
flow  rate  demands  of  5  and  10  SLPM  (standard  liter  per 
minute).  The  results  in  Fig.  2  show  that  the  dynamic  pres¬ 
sure  behavior  is  predicted  quite  well  by  the  RC  form  of  the 
model,  especially  during  the  periods  where  the  pressure  re¬ 
mains  relatively  constant.  The  minor  discrepancies  observed 
just  before  the  pressure  begins  to  level  off  are  most  likely 


associated  with  the  one-dimensional  model  not  accounting 
for  the  convection  of  energy  in  the  axial  direction,  as  noted 
earlier.  This  becomes  more  apparent  in  the  prediction  of  the 
temperature  histories  shown  in  Fig.  3.  The  RC  form  of  the 
model  does  quite  well  in  predicting  the  temperature  history 
for  the  5  SLPM  hydrogen  flow  rate  demand,  but  begins  to  de¬ 
viate  a  bit  more  at  the  higher  hydrogen  flow  rate  demand  of  10 
SLPM,  but  still  less  than  1  K.  It  can  therefore  be  stated  that, 
for  hydrogen  flow  rate  demands  on  the  order  of  10  SLPM  or 
less,  where  the  heat  effects  are  not  so  drastic,  the  RC  form  of 
the  model  is  quite  sufficient  for  use  to  study  the  performance 
of  dynamically  coupled  systems. 


3.  Proton  exchange  membrane  (PEM)  fuel  cell  model 


The  main  power  source  of  the  integrated,  hydrogen-fueled 
fuel  cell  system  is  a  25 -cell  PEM  fuel  cell  stack  with  the  cells 
organized  in  series.  The  voltage  of  each  cell  varies  with  the 
current,  temperature  and  pressure  according  to  the  empirical 
equation  shown  below  [18]. 


^stack  —  N 


Eo  ~  (bo TftQ  +  b\T^c  +  hi)  log 


A  cell 


-  (7*0 TpC  +  n^FC  +  r2))~ - (^O^FC  +  ml^FC 


A  cell 


+  m2)  exp  ( (fto^FC  +^1^FC  +^2) 


Acell  /  _ 


(22) 


where  Eq  is  the  standard  potential  of  H2-O2  reaction,  I  the 
current  from  the  fuel  cell,  £stack  the  voltage  of  the  fuel  cell 
stack,  N  the  number  of  cells,  Aceii  the  area  of  each  cell,  Tpc 
the  bulk  temperature  of  the  fuel  cell  stack,  and  Ppc  is  the 
hydrogen  pressure  in  the  fuel  cell.  Constants  bu  77,  mu  and 
rii  (i  =  0,  1,2)  are  coefficients  determined  empirically  from 
experimental  data.  This  fuel  cell  model  was  also  validated 
against  experimental  data  elsewhere  [18]  (also  refer  to  the 
online  document  [11]). 

The  relationship  between  the  hydrogen  molar  flow  rate  h 
and  the  current  I  generated  by  the  fuel  cell  is  given  by 


2  hF 

I  =  - 

N 


(23) 


where  F  is  the  Faraday  constant.  The  thermal  power  W± 
generated  by  the  fuel  cell  is  given  by  the  following  formula: 


Wth  =  hAHc  -  WQ 


(24) 


where  A Hc  is  the  heat  of  combustion  of  hydrogen  and  WQ  the 
electrical  power  generated  by  the  fuel  cell.  Also,  the  thermal 
power  released  from  the  fuel  cell  stack  is  divided  into  two 
parts:  some  is  transferred  to  the  metal  hydride  bed,  and  the 
rest  is  transferred  to  the  ambient  environment.  This  situation 
is  denoted  by: 


Fig.  3.  Model  predictions  of  experimental  temperature  histories  [10]  in  the 
metal  hydride  bed  for  hydrogen  flow  rate  demands  of  5  and  10  SLPM. 


With  =  Q  +  Wamb  +  Cpc^FC 


dTpc 
d  t 


(25) 
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where  Cpc  is  the  average  heat  capacity  of  the  fuel  cell  stack, 
rape  is  the  mass  of  the  fuel  cell  stack,  and  Wa mb  is  the  thermal 
power  transferred  between  the  fuel  cell  heat  exchanger  and 
the  ambient  environment  according  to 


fkamb  —  ^ex^ex(^FC  ^amb)  (26) 

where  hQX  is  the  overall  heat  transfer  coefficient  of  the  fuel 
cell  heat  exchanger  and  Aex  is  the  effective  area  of  the  heat 
exchanger.  The  fuel  cell  stack  voltage  and  current  are  related 
to  the  electrical  power  generated  by  the  fuel  cell  by 

We  =  £stack/  (27) 


4.  System  level  simulation  in  the  VT  environment 

System  level  simulations  of  the  hydrogen-fueled,  fuel  cell 
system  were  carried  out  with  and  without  the  fuel  cell  be¬ 
ing  thermally  coupled  to  the  metal  hydride  bed.  In  both  sit¬ 
uations  the  metal  hydride  bed  is  connected  to  the  fuel  cell 
through  a  pressure  regulator  to  supply  hydrogen  at  1 .0  atm 
constant  pressure.  Air  is  also  fed  to  the  fuel  cell  at  1 .0  atm 
constant  pressure.  For  the  uncoupled  simulations,  both  the 
fuel  cell  and  the  metal  hydride  bed  are  configured  to  indepen¬ 
dently  exchange  heat  with  the  ambient  environment  (i.e.,  the 
temperature  source).  This  integrated,  but  thermally  uncou¬ 
pled  hydrogen  storage  and  fuel  cell  energy  system  is  shown 
schematically  in  the  VTB  environment  in  Fig.  4a.  The  water 
circulating  through  the  metal  hydride  bed  and  the  fuel  cell 
heat  exchangers  (labeled  1  and  2,  respectively)  in  this  case 
is  held  constant  at  296  K,  the  ambient  temperature.  For  the 
coupled  simulations,  the  metal  hydride  bed  is  configured  to 
exchange  energy  with  the  fuel  cell  stack  through  a  heat  ex¬ 
changer,  since  the  fuel  cell  produces  excess  (waste)  heat  dur¬ 
ing  operation.  Fig.  4b  shows  the  VTB  environment  schematic 
of  this  integrated  and  now  thermally  coupled  hydrogen  stor¬ 
age  and  fuel  cell  energy  system.  Except  for  the  heat  transfer 
between  the  metal  hydride  bed  and  the  fuel  cell  being  cou¬ 
pled,  the  rest  of  the  system  is  identical  to  that  shown  in  Fig.  3. 
However,  in  this  case,  the  temperature  of  the  recirculating 
water  is  assumed  to  be  equal  to  the  bulk  temperature  of  the 
fuel  cell.  The  values  of  all  the  parameters  used  in  the  model 
equations  are  given  in  Tables  1  and  2. 

For  the  thermally  coupled  case,  1 1  equations,  namely  Eqs. 
(4),  (14)— (17),  and  (22)-(27),  were  solved  simultaneously  for 
11  unknowns,  namely,  Es tack,  l  n,  P ,  q ,  q\  g,  T,  7pC,  ITamb 
and  With*  Also,  for  this  coupled  case,  the  temperature  of  the 
metal  hydride  bed  heat  exchanger  medium  Tm  in  Eq.  (17) 
becomes  Tpc-  For  the  uncoupled  case,  the  same  set  of  equa¬ 
tions  and  unknowns  are  solved  simultaneously;  however,  Eq. 
(25)  is  no  longer  valid  because  g  does  not  exchange  energy 
with  the  fuel  cell  but  instead  with  the  ambient  environment. 
In  this  case,  Eq.  (25)  simply  becomes  W±  =  ITamb  and  the 
temperature  of  the  metal  hydride  bed  heat  exchanger  medium 
Tm  in  Eq.  (17)  becomes  Tamt,. 


Temperature  Source 


Metal  Hydride  Bed  Pressure  Regulator 

i  i  i  I  PMF  Fuel  Cell 


. .  ”  Pressure  Reference 

Heat  Exchanger  1 


Pulse  Power  Load 
"  Electrical  Ground 


Temperature  Source 


Fig.  4.  Schematics  of  thermally  uncoupled  and  thermally  coupled  hydro¬ 
gen  storage  and  fuel  cell  systems  in  the  VTB  environment:  (a)  thermally 
uncoupled  system  and  (b)  thermally-coupled  system. 


In  both  the  uncoupled  and  coupled  simulations,  a  pulse 
power  load  (ITe)  draws  200  W  for  0.25  h  and  then  100  W  for 
0.25  h  alternately  from  the  fuel  cell  for  4  h,  as  shown  in  Fig.  5. 
The  corresponding  voltage  and  temperature  of  the  fuel  cell 
stack  are  given  in  Figs.  6  and  7,  respectively.  Fig.  8  shows  the 
thermal  power  released  from  the  fuel  cell  during  operation  of 


Fig.  5.  Pulsed  power  drawn  from  the  fuel  cell  stack,  alternating  between 
200  W  and  100  W  every  0.5  h  for  4  h  of  operation,  whether  thermally  coupled 
to  the  metal  hydride  bed  or  not. 
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Table  1 


Model  parameters  used  in  the  fuel  cell  and  metal  hydride  bed  simulations 
[10] 


Parameter 

Value 

Aex  (cm2) 

2.00E+2 

A  cell  (cm2) 

2.00E+2 

bio  (Pa- 1/2 ) 

2.36E-6 

b2o  (Pa- 1/2 ) 

7.87E-6 

CFC*  (JK_1  kg-1) 

8.13E+2 

Cph2  (J  mol- 1  K"1) 

2.88E+1 

Cps  (Jmol-'K-1) 

4.20E+2 

Do 

2.40E+1 

Eo  (V) 

1.00E+0 

F  (C  mol-1) 

9.65E+4 

Ktd  (W  m“2  K”1) 

7.57E+1 

hex  (W  m-2  K-1) 

1.70E+2 

kq  (s'1) 

8.14E-3 

E  (m) 

1.52E+0 

(kg) 

1.72E+0 

MMh  (kg  mol-1) 

7.29E-2 

Po  (Pa) 

5.05E+6 

Pfc  (Pa) 

1.01E+5 

9WX  (mol  mol- 1 ) 

9.84E-1 

qs  (mol  mol-1) 

1.25E-1 

R  (J  mol-1  K-1) 

8.31E+0 

Ri  (m) 

2.00E-2 

R0  (m) 

4.50E-2 

Tamb  (K) 

2.96E+2 

Pc  (K) 

2.96E+2 

To  (K) 

2.96E+2 

£ 

5.28E-1 

Ps  (kgm-3) 

7.61E+3 

N 

2.50E+1 

AH\  (Jmol-1) 

—  1.38E+4 

AH2  (J  mol- 1 ) 

-206E+4 

AHS  (J  mol-1) 

— 2.77E+4 

A Hc  (J  mol-1) 

— 2.42E+5 

*  Based  on  [19,20]. 


Fig.  6.  Voltage  of  the  fuel  cell  stack  during  pulsed  operation  of  the  fuel  cell 
while  thermally  coupled  and  uncoupled  to  the  metal  hydride  bed. 


Fig.  7.  Temperature  of  the  fuel  cell  stack  during  pulsed  operation  of  the  fuel 
cell  while  thermally  coupled  and  uncoupled  to  the  metal  hydride  bed. 

the  fuel  cell  over  the  4  h  for  both  uncoupled  and  coupled  op¬ 
erations.  The  hydrogen  flow  rate  demand  and  the  cumulative 
percentage  of  the  hydrogen  in  the  metal  hydride  bed  con¬ 
sumed  over  this  4h  period  are  shown  in  Fig.  9.  Figs.  10-12, 
respectively,  display  the  corresponding  hydrogen  loading  in 
the  solid  metal  hydride  phase,  the  metal  hydride  bed  pressure, 
and  the  metal  hydride  bed  temperature  histories  for  both  the 
uncoupled  and  coupled  systems.  Finally,  Fig.  13  shows  the 
thermal  power  transferred  to  the  metal  hydride  bed  from  the 
fuel  cell  over  the  4h  for  both  the  uncoupled  and  coupled 
operations. 

During  the  4  h  period  of  pulsed  operation  of  the  fuel  cell, 
with  a  load  alternately  drawing  200  W  and  100W  every 
0.25  h,  the  fuel  cell  operates  at  nearly  a  constant  voltage  of 
22.5  zb  0.1  V  (Fig.  6).  The  0.1  V  periodic  swing  is  due  to  the 
high  and  low  power  levels,  with  the  higher  voltage  naturally 
corresponding  to  the  lower  power  level.  The  uncoupled  sys¬ 
tem  also  operates  at  a  marginally  higher  voltage,  most  likely 
due  to  the  effect  of  thermal  coupling  on  the  fuel  cell  temper¬ 
ature  during  pulsed  operation. 

Fig.  7  shows  that  during  this  4  h  period,  the  fuel  cell  tem¬ 
perature,  initially  starting  from  296  K,  exhibits  essentially  a 

Table  2 


Values  of  the  polynomial  coefficients  used  in  the  PEM  fuel  cell  model  [18] 


Coefficient 

Value 

d0  (mV  K-2  cell-1) 

3.09E-3 

di  (mV  K-1  cell-1) 

— 2.05E+0 

d2  (mV  cell-1) 

4.60E+2 

ro  (£2  cm2  K-2  cell-1) 

1.52E-5 

r\  (£2  cm2  K-1  cell-1) 

—  1.00E— 2 

r2  (£2  cm2  cell-1) 

2.50E+0 

mo  (mV  Pa-2  cell-1) 

1 .3  IE —  1 1 

m\  (mV  Pa-1  cell-1) 

-7.33E-6 

m2  (mV  cell-1) 

1.05E+0 

no  (cm2  mA-1  Pa-2) 

— 2.84E— 13 

n\  (cm2  mA-1  Pa-1) 

1.16E-7 

n2  (cm2  mA- 1 ) 

-6.39E-3 
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Time  (hour) 

Fig.  8 .  Thermal  power  released  from  the  fuel  cell  to  the  ambient  environment 
during  uncoupled  operation,  and  to  the  ambient  environment  and  the  metal 
hydride  bed  during  coupled  operation. 

saw  tooth  pattern  and  oscillates  between  304  K  and  300  K 
corresponding  to  the  high  and  low  power  periods,  respec¬ 
tively.  It  is  noteworthy  that  a  slightly  higher  temperature  is 
observed  for  the  uncoupled  operation  especially  during  the 
high  power  periods.  This  suggests  that  more  energy  is  be¬ 
ing  removed  from  the  fuel  cell  when  the  metal  hydride  bed 
is  thermally  coupled  to  the  fuel  cell,  a  very  positive  but  in 
this  case  a  very  minor  effect.  This  result  is  further  substan¬ 
tiated  in  Fig.  8,  which  shows  that  during  coupled  operation, 
slightly  more  energy  is  transferred  from  the  fuel  cell  during 
high  power  operation.  Fig.  8  also  shows  that  during  the  4  h 
period  of  pulsed  operation,  with  the  load  varying  between 
200  W  and  100  W,  approximately  80  W  and  40  W  of  power 
are  respectively  converted  into  waste  heat.  This  means  that 
60%  of  the  energy  obtained  from  the  hydrogen  is  converted 
into  useful  electricity.  The  other  40%  of  the  energy  is  con¬ 
verted  into  waste  heat,  with  no  other  energy  losses  accounted 


Fig.  9.  Hydrogen  flow  rate  demand  and  cumulative  percentage  of  initial 
hydrogen  in  the  metal  hydride  bed  consumed  over  the  4  h  period  of  pulsed 
operation  of  the  fuel  cell  whether  coupled  to  the  metal  hydride  bed  or  not. 


Fig.  10.  Hydrogen  loading  in  the  solid  metal  hydride  phase  in  the  metal 
hydride  bed  during  pulsed  operation  of  and  while  thermally  coupled  and 
uncoupled  to  the  fuel  cell. 

for  in  this  system  simulation.  This  brings  up  one  of  the  main 
questions  addressed  in  this  work:  what  is  the  effect  of  trans¬ 
ferring  some  of  this  waste  energy  from  the  fuel  cell  to  the 
metal  hydride  bed. 

Fig.  9  shows  that  during  the  4  h  of  pulsed  operation  the  fuel 
cell  consumes  a  little  more  than  50%  of  the  total  hydrogen  in 
the  bed,  with  the  hydrogen  flowing  at  constant  but  alternating 
rates  of  1.2  and  0.6  SLPM,  respectively  and  proportionally, 
corresponding  to  the  high  and  low  power  levels.  However, 
the  hydrogen  consumption  is  the  same  whether  the  metal  hy¬ 
dride  bed  is  thermally  coupled  to  the  fuel  cell  or  not,  which 
is  clearly  reflected  in  the  change  in  the  hydrogen  loading  in 
the  metal  hydride  phase.  Fig.  10  shows  that  it  decreases  from 
about  0.81  mol  mol-1  to  0.43  mol  mol-1 .  Hence,  a  little  less 
than  half  of  the  hydrogen  in  the  solid  (metal  hydride)  phase 
is  consumed  during  4h  of  operation;  accounting  also  for  the 


xio5 


Fig.  11.  Gas  phase  pressure  inside  the  metal  hydride  bed  during  pulsed 
operation  of  and  while  thermally  coupled  and  uncoupled  to  the  fuel  cell. 
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Fig.  12.  Average  bulk  temperature  of  the  metal  hydride  bed  during  pulsed 
operation  of  and  while  thermally  coupled  and  uncoupled  to  the  fuel  cell. 

hydrogen  in  the  gas  phase  brings  the  total  hydrogen  con¬ 
sumed  to  a  little  more  than  50%  as  stated  above  and  shown 
in  Fig.  9.  So,  as  may  be  expected,  thermally  coupling  the 
system  does  not  affect  the  hydrogen  flow  rate,  which  is  dic¬ 
tated  by  the  power  requirements  of  the  load  imparted  to  the 
fuel  cell.  It  also  does  not  affect  the  hydrogen  loading  in  the 
metal  hydride  phase,  a  result  that  is  not  so  obvious  but  is 
easily  rationalized  because  this  is  the  source  of  the  hydrogen. 
What  thermally  coupling  does  effect  is  the  temperature  of  the 
metal  hydride  bed,  which  in  turn  effects  the  pressure  of  the 
bed.  These  effects  are  clearly  observed  in  Figs.  11  and  12. 
The  different  curves  in  these  figures  reflect  the  differences 
observed  when  the  metal  hydride  bed  exchanges  energy  with 
the  ambient  environment  (uncoupled)  compared  to  when  it 
exchanges  energy  with  the  fuel  cell  (coupled),  as  shown  in 
Fig.  13. 


Time  (hour) 


Fig.  13.  Thermal  power  transferred  to  the  metal  hydride  bed  from  the  am¬ 
bient  environment  during  uncoupled  operation  and  from  the  fuel  cell  during 
coupled  operation. 


For  the  uncoupled  case,  the  pressure  inside  the  bed  de¬ 
creases  from  about  5.0  x  105  Pa  to  about  3.3  x  105  Pa.  How¬ 
ever,  in  the  coupled  case,  the  pressure  inside  the  bed  decreases 
much  less  from  about  5.0  x  105  Pa  to  only  about  4. 1  x  105  Pa 
over  the  same  4  h  of  operation.  This  is  due  to  the  more  en¬ 
ergy,  on  average,  being  transferred  to  the  metal  hydride  bed 
when  the  bed  is  thermally  coupled  with  the  fuel  cell  (about 
3  W,  net)  compared  to  when  it  is  thermally  coupled  with  the 
ambient  environment  (about  1  W  net),  as  shown  in  Fig.  13. 
This  small  amount  of  extra  energy  causes  the  temperature 
and  hence  the  pressure  of  the  metal  hydride  bed  to  both  in¬ 
crease,  with  the  temperature  cycling  between  297.6  K  and 
299.2  K.  Although  these  temperatures  are  always  above  the 
ambient  temperature  of  296  K,  this  is  not  the  temperature  of 
interest  for  the  coupled  case;  the  temperature  of  interest  is  the 
fuel  cell  temperature,  which  varies  between  300  K  and  304  K 
(Fig.  7).  Hence,  the  metal  hydride  bed  is  always  operating 
below  the  fuel  cell  temperature,  which  means  that  energy  is 
always  being  transferred  to  the  bed  (Fig.  13),  a  positive  ef¬ 
fect.  In  contrast,  for  the  uncoupled  case,  the  metal  hydride 
bed  temperature  is  always  below  the  ambient  temperature, 
cycling  between  291.3  K  and  293.5  K.  In  this  case,  the  ambi¬ 
ent  temperature  is  the  temperature  of  interest,  and  again  with 
the  bed  always  operating  below  the  ambient  temperature,  en¬ 
ergy  is  being  transferred  to  the  bed  (Fig.  13).  However,  the 
higher  operating  temperature  of  the  thermally  coupled  metal 
hydride  bed  translates  into  a  greater  time  during  which  the 
fuel  cell  may  continue  to  produce  200  W  and  100  W  of  useful 
power.  This  is  due  to  the  pressure  staying  above  the  lowest 
delivery  pressure  of  the  fuel  cell  for  a  longer  period  of  time, 
in  this  case  1  atm.  Of  course,  eventually  the  metal  hydride 
bed  will  begin  to  run  out  of  hydrogen  and  the  pressure  will 
decrease  below  this  minimum  operating  pressure  no  matter 
how  high  the  temperature. 

It  is  interesting  that  the  temperatures  of  the  metal  hydride 
bed  for  the  coupled  and  uncoupled  cases  both  exhibit  saw¬ 
tooth  patterns  but  oscillate  1 80°  out  of  phase  with  each  other, 
i.e.,  when  the  temperature  is  increasing  in  one  case  it  is  de¬ 
creasing  in  the  other  case  and  vice  versa.  Not  surprisingly, 
the  thermal  energy  imparted  to  these  beds  also  shows  these 
opposite  behavior,  but  with  drastically  different  patterns.  For 
the  uncoupled  case,  the  temperature  of  the  metal  hydride  bed 
initially  decreases  from  296  K  to  291  K  during  the  first  high 
power  cycle;  thereafter,  it  exhibits  2  K  swings  in  tempera¬ 
ture  between  293  K  and  291  K,  always  being  below  the  am¬ 
bient  temperature  of  296  K.  Hence,  energy  is  always  being 
transferred  to  the  bed,  a  favorable  effect  as  explained  below. 
The  periodic  decreases  in  temperature  occur  when  the  pulsed 
power  source  draws  200  W  from  the  fuel  cell;  in  contrast, 
when  only  100  W  of  power  are  drawn,  the  temperature  in¬ 
creases.  The  opposite  behavior  is  exhibited  by  the  thermal 
power  transferred  to  the  metal  hydride  bed  from  the  ambient 
environment.  Initially,  about  2.7  W  are  transferred  at  most  to 
the  metal  hydride  bed  and  then  the  thermal  power  oscillates 
between  2.7  W  and  1.5  W  when  the  temperature  of  the  metal 
hydride  bed  is  lowest  and  highest,  respectively.  This  coupled 
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and  1 80°  out  of  phase  behavior  between  the  temperature  of 
and  the  thermal  energy  transferred  to  the  bed  is  explained 
by  the  metal  hydride  bed  cooling  when  it  releases  hydrogen 
[12].  The  greater  this  decrease  in  temperature,  the  greater 
the  driving  force  for  heat  transfer  between  the  bed  and  the 
ambient  environment,  the  greater  the  power  input  to  the  bed. 
So  when  the  bed  periodically  heats  back  up  to  293  K,  the 
driving  force  for  heat  transfer  decreases,  which  causes  the 
power  input  to  decrease.  Nevertheless,  the  energy  transferred 
to  the  bed  when  only  100  W  are  being  drawn  from  the  fuel 
cell  is  apparently  enough  to  overcome  the  endothermic  pro¬ 
cess  associated  with  the  release  of  hydrogen  at  the  lower  flow 
rate  demand  and  cause  it  to  heat  back  up  to  293  K.  However, 
when  the  power  source  draws  200  W  and  the  hydrogen  flow 
rate  demand  correspondingly  increases  to  the  higher  value, 
although  the  energy  transferred  to  the  bed  increases,  it  is  not 
enough  to  overcome  the  endothermic  process  associated  with 
the  release  of  this  increased  amount  of  hydrogen.  Hence,  the 
bed  cools  during  this  period  to  about  29 1  K.  This  periodic 
behavior  continues  over  the  4  h  of  operation  of  the  fuel  cell. 

Similar  temperature  change  and  corresponding  driving 
force  arguments  are  made  to  describe  the  behavior  of  the  cou¬ 
pled  system.  However,  in  this  case,  the  3. 7-2. 7  W  of  power 
imparted  to  the  bed  during  high  power  operation  is  apparently 
enough  to  offset  the  endothermic  cooling  effect  mentioned 
above  because  the  temperature  of  the  bed  actually  increases 
during  this  part  of  the  cycle.  This  necessarily  causes  the  ther¬ 
mal  power  transferred  to  the  bed  to  decrease  during  the  same 
period,  because  although  the  fuel  cell  temperature  is  at  its  high 
value,  it  essentially  remains  constant  during  this  period.  At 
the  onset  of  low  power  operation,  the  temperature  of  the  fuel 
cell  decreases  instantaneously  to  its  low  value  (Fig.  7),  which 
lowers  the  thermal  power  input  to  the  metal  hydride  bed  to 
about  0.5  W  due  to  the  small  driving  force  for  heat  transfer. 
This  allows  the  endothermic  effect  to  begin  to  dominate  even 
though  the  hydrogen  demand  is  low  and  the  temperature  of 
the  metal  hydride  bed  begins  to  decrease.  Nevertheless,  the 
thermal  power  input  to  the  bed  increases  during  this  period  to 
about  1.5  W  because  the  driving  force  for  heat  transfer  now 
increases  due  to  the  decrease  in  the  metal  hydride  bed  tem¬ 
perature  with  the  fuel  cell  temperature  essentially  remaining 
constant.  This  periodic  behavior  continues  over  the  4h  of 
operation  of  the  fuel  cell;  however,  the  periodic  behavior  is 
clearly  and  subtly  different  when  compared  to  the  uncoupled 
case. 


5.  Conclusions 

The  overall  objective  of  this  paper  was  to  study  the  be¬ 
havior  of  a  thermally  coupled  hydrogen  storage  and  fuel  cell 
system.  An  underlying  and  more  specific  objective  was  to 
determine  the  effect  of  transferring  waste  (thermal)  energy 
from  the  fuel  cell  to  the  metal  hydride  bed  in  a  truly  dynamic 
and  coupled  fashion.  These  objectives  were  met  by  formu¬ 
lating  dynamic  models  of  a  metal-hydride  hydrogen  storage 


system  and  a  PEM  fuel  cell  stack  for  easy  implementation 
in  the  virtual  test  bed  (VTB)  environment,  a  simulation  plat¬ 
form  designed  for  carrying  out  systems  level  modeling  like 
the  integrated  fuel  cell  system. 

An  inherently  nonlinear  model  that  realistically  accounts 
for  heat  and  mass  transfer  resistances  associated  with  the 
discharge  of  hydrogen  from  a  metal  hydride  bed  was  refor¬ 
mulated  in  a  resistive  companion  (RC)  form  and  success¬ 
fully  validated  against  experimental  discharge  data.  The  RC 
form  of  the  metal  hydride  bed  model  was  then  coupled  with 
a  relatively  simple,  experimentally  validated,  PEM  fuel  cell 
model.  These  coupled  models  were  used  to  carry  out  dynamic 
simulations  of  the  integrated  fuel  cell  system  in  the  VTB  en¬ 
vironment.  Two  specific  situations  were  studied:  in  the  first 
scenario,  the  fuel  cell  and  metal  hydride  bed  were  thermally 
coupled;  in  the  second  scenario,  they  were  thermally  uncou¬ 
pled. 

The  results  from  these  simulations  clearly  revealed  unique 
and  subtle  behavior  associated  with  thermally  coupled  sys¬ 
tems  that  could  not  be  easily  gleaned  from  simulating  each 
device  alone  as  revealed  from  studying  the  uncoupled  case. 
For  example,  the  results  exemplified  the  importance  of  pro¬ 
viding  energy  to  the  metal  hydride  bed  to  facilitate  the  re¬ 
moval  of  hydrogen  during  a  discharge  process.  Although  a 
substantial  amount  of  energy  in  the  form  of  waste  heat  tends 
to  be  available  from  a  PEM  fuel  cell  to  increase  the  temper¬ 
ature  of  the  metal  hydride  bed  during  discharge,  this  work 
showed  that  in  some  cases,  very  little  energy  may  be  needed 
to  make  a  substantial  difference  in  the  performance  of  the 
metal  hydride  bed.  A  relatively  small  increase  in  the  operat¬ 
ing  temperature  of  the  metal  hydride  bed,  which  occurred  for 
the  thermally  coupled  scenario,  had  a  significant  and  positive 
effect  on  increasing  the  pressure  in  the  metal  hydride  bed. 
This  sustained  higher  pressure  necessarily  improves  the  uti¬ 
lization  of  stored  hydrogen  because  the  pressure  of  the  bed 
is  able  to  stay  above  the  lowest  operating  pressure  of  the  fuel 
cell  for  a  longer  period  of  time.  This,  in  turn,  translates  into  a 
longer  operating  time  for  the  fuel  cell.  A  detailed  parametric 
study  is  forthcoming. 
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